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ABSTRACT 


An investigati:. vas performed to evaluate the accuracy 
of numerical mode. ~ of near-fields for antennas on or near 
surfaces using =. Numerical Electromagnetic Code (NEC). 


Average power gain and input impedance were calculated, for 
two models. The first, a dipole antenna located inside a 
cube-shaped wire-grid box of 1 meter sides was evaluated for 
a wide range of frequencies in free space. The second, a 
monopole antenna mounted on the top of two cube shaped boxes 
(a wire grid and a surface patch box) of .1 meter sides over 
a perfectly conducting ground was evaluated from 1 to 1.4 
GHZ. The monopole was positioned at the center, at an edge, 
and at the corner of the box top. For the dipole in the 
box, the average gain and input impedance are presented and 
evaluated. For the monopole on the box, calculated results 
are compared to measured admittance values. The NEC wire- 
grid model results compare closely to the measurements, but 
for surface-patch models, only one position of the monopole 
yielded satisfactory correlation to measurements. Recommen- 
dations for improved numerical modeling performance are 


made. 
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I. INTRODUCTION 


A. GENERAL BACKGROUND 

Because shipboard operations are carried out within 
fixed (small) distances from transmitting antennas, the Navy 
has a unique and long-standing operational problem: the 
radiation from these antennas can be dangerous to personnel, 
ordnance, fuel, and electronic equipment. This is called 
RADIATION HAZARD (RADHAZ). 

Accordingly, the Navy has been pursuing the study of 
near fields of antennas for a number of years. However, the 
near field structure is complex, and previous theoretical 
analysis has been practical only for simple antennas in 
uncomplicated geometrical settings. With the advent of the 
modern high-speed computer, approximate solution techniques 
such as the method of moments used by the Numerical Electro- 


magnetic Code have become practical.[{Ref. 1] 


B. NEED FOR STUDY 

With the added emphasis on near field (due to RADHAZ) by 
EPA and USN demands, we must be more accurate in our near 
field predictions to assure that personnel and equipment 
can continue to function = safely in modern shipboard 


RF environments. 


C. SCOPE AND LIMITATIONS 

Two theoretical models of cube shaped boxes (similar to 
those found on = structures on shipboard systems) are 
examined. These boxes are modeled numerically and the 
parameters of one of the models is compared to measured 
values from Reference 2. Modeling techniques are varied to 
assess their impact on accuracy. Guidelines are developed 
for maximum accuracy. Also, it is known that the Numerical 
Electromagnetic Code has been found to suffer loss of 
accuracy in VLF applications involving electrically small 
antennas [Refs. 3 and 4]. This may affect the modeling 


efforts of Ehis thesics 


II. NEAR FIELD THEORY 


A. EXTERIOR FIELDS OF RADIATING ANTENNAS 

Figure 2.1 describes the regions into which the exterior 
fields of a radiating antenna are commonly divided. The 
antenna radiates into free space as a linear system with the 
single-frequency time sare nionee of e(-jut), 

The FAR-FIELD region extends to infinity and is that 
region of space where the radial dependence of electric and 
magnetic fields varies approximately as e(-JKY) yr. The 
inner radius of the far field can be estimated from the 
general free-space integral for the vector potential and is 
usually set at DN +(2D2/)) . (The added »& covers the pos- 
sibility of the maximum dimension D of the antenna being 
smaller than a wavelength. In other words, the RAYLEIGH 
distance (z)t should actually be measured from the outer 
boundary of the reactive near field of the antenna.) For 
the main beam direction, the RAYLEIGH distance can sometimes 
be reduced. 

The free-space region from the surface of the antenna to 
the far field is labeled as the NEAR-FIELD region. It is 


divided into two subregions, the reactive and the 


1 The Rayleigh distance is given by Zg = mD?/) 





RADIATING FAR FIELD 
NEAR FIELD Fiat) aia 
/| 
ea 
| 
lL [Fresnel | FRAUNHOFER 
FIELD | ! | | 
DISTRIBUTION ———> | | 
i : 
ow . 
2 | fee tees 
0 A So! OD Pisione RANGE fr 
Te 


Figure 2.1 Exterior fields of radiating antennas [Ref.5]. 


radiating near fields. The reactive near-field region is 
commonly thought to extend about A/27 from the surface of 
the antenna, although experience with near-field measure- 
ments indicates that a distance of a wavelength (A) or so 
would be a more reasonable outer boundary to the reactive 
near field. 

The reactive near field can be defined in terms of 
planar, cylindrical, or spherical modes. Unfortunately, the 
reactive fields of spherical (or cylindrical) multipoles are 
not identical to the plane-wave evanescent fields of 
multipoles. Thus, a simpler, physically appealing and less 


ambiguous method of defining the reactive region of antennas 


relies directly on Poynting's theorem and the vector 
potential. It can be shown that the contribution to the 
reactive part of the input impedance of an antenna from the 
fields outside a surface surrounding the antenna is 
proportional to the imaginary part of the complex Poynting 
vector integrated over the surface. Thus, wherever the 
phase of the electric and the magnetic field vectors are 
near quadrature, the Poynting vector will contribute mainly 
to the reactive part of the input impedance. Taking the 
curl of the vector potential integral once to get the 
Magnetic field and twice to get the electric field, shows 
that the phase of the electric and magnetic fields may be 
near quadrature in regions within approximately a wavelength 
(A) of the antenna. Consequently, the region within about a 
wavelength of the physical antenna is referred to as the 
REACTIVE NEAR FIELD. 

Beyond a distance of about a wavelength, the electric 
and magnetic fields tend to propagate predominantly in 
phase, but, of course, do not exhibit e (-JKYr) yy dependence 
until the far field is reached. This propagating region 
between the reactive near field and the far field is called 
the RADIATING NEAR FIELD. 

Finally, the optical terms "Fraunhofer region" and 
"Fresnel region" are sometimes used to characterize the 
fields of antennas. The term FRAUNHOFER REGION can be used 


synonymously with the far-field region or to refer to the 


focal region of an antenna focused at a finite distance. 
The FRESNEL REGION, which extends from about (D/2)1/3 pD/2 
+ A to the far field, is the region up to the far field in 
which a quadratic phase approximation can be used in the 
vector potential integral. The Fresnel region is a sub- 


region of the radiating near-field region.[Ref. 5:pp. 33-34] 


B: DIFFICULTIES WITH NEAR FIELD CALCULATIONS 

The NEAR FIELD is more difficult to calculate than the 
Far field. When calculating the near zone radiation, the 
antenna no longer appears as an infinitesimal point (as is 
the assumption for far field), but rather we are conscious 
that the energy coming toward us comes from separated 
locations. Also the terms of the field expressions with 
powers of 1/r (r is the distance from the origin of the 
antenna to the field point) greater than one are here more 
important than the 1/r-dependent terms. The complex 
Poynting vector L(E X H*) will then contain terms with a 
power of (1/r) in addition to the radiation field term. It 
results that these terms are purely imaginary, indicating 
reactive power and energy oscillating in and out. The near 
field is restricted to charge and current density and far 
field only to current. All these difficulties make the near 
field calculation more imprecise. 

In the past, accuracy of near-field calculations was not 
a problem; there was no practical application to which the 
energy could be applied. However, in recent years there has 
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been an escalating interest in the close-in problen, 
especially in the Navy where, in general, antennas mounted 
on a ship will couple into metallic parts of the ship. This 
means that ecrena RF currents can be excited on closely 
coupled conducting surfaces, which in turn re-radiate. The 
total near field is the (vector) sum of the fields radiated 
by the antenna, the primary antenna, and the nearby 
conducting BRE eCe: Re-radiation is enhanced whenever the 
secondary scatterers are similar in size to the prymary 
antennas or whenever resonant length conducting paths and 
spacing occur. Because of this the near field produces 
RADIATION HAZARD (RADHAZ) that can affect personnel, 


ordnance, and fuel on ships. 


LIT. NUMERICAL MODELS OF SURFAGES TY PRESS TRUGCRURES 


A. NUMERICAL MODELING OF SURFACES FOR NEAR-FIELD 
PREDICTION 


There are two principal ways in which the method of 
moments can be used to model antennas together with the 
Surrounding structures. The simplest way to model a solid 
surface is with a grid of wires, the so-called wire-grid 
model. The other approach is to divide the surface into 
patches or cells, each having a continuous metallic surface. 
This is the so-called surface-patch model. 

The successful substitution of a wire grid or surface 
patch for a continuous metallic surface depends upon the 
fact that as the model (the wire grid or surface patch) 
becomes smaller relative to a wavelength, it supports a 
current distribution which approximates that of the 
continuous surface. The current is only an approximation to 
the actual current, however, and as such it can be expected 
to accurately predict the far fields but not necessarily the 
near fields. This is due to the fact that the wire-grid 
model supports an evanescent reactive field on both sides of 
its surface. An actual continuous conducting surface is not 
capable of supporting such a field. For surface patches the 
current is assumed to be continuous over a patch which is 
only a simple approximation to the current and allows no 
smooth variation over the body being modeled. 
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The version of the Numerical Electromagnetic Code that 
will be used in this work is NEC-3, although it has been 
reported by Burke [Refs. 3 and 4] that there are some 
inaccuracies in NEC. [In the first model, the wire grid with 
interior dipole, there are some frequencies in which the 
wavelengths, are more than ten-times larger than the side 
dimension of the cell segments, and this may affect results. 
To increase confidence in the accuracy of the results, two 
antenna parameters are calculated: the average power gain 
and the input impedance. The average power gain provides a 
good check of the solution accuracy if it is close to one 
for free space or close to two for a perfectly conducting 
ground plane. For our purposes we will consider average 
power gain within +10% of theoretical to be adequate for 
engineering purpose. Correct average power gain does not 
insure accurate near-field prediction. Therefore, input 
impedance, which is in general complex, is also checked. 
The real part of the impedance relates to radiated far-field 
power while the imaginary portion relates to the balance of 
reactive energy in the near field. 

In the first model analyzed, the wire grid with a 
dipole inside, there are no measured values for comparison. 
In the second model, the wire-grid box with the monopole, 
the measured admittance is available. The next section 


describes these two models. 


B. WIRE GRID BOX WITH INTERIOR DIPOLE 

A WIRE-GRID box with an interior dipole was modeled to 
obtain numerical values of average power gain and input 
impedance. The model constructed was a 2 meter per side 
cube-shaped box, with a dipole inside, 1 Meter long. The 
grid was .4 by .4 meter cells, with each cell side divided 
into two segments (see Figure 3.1). (Initial wire grid box 
models employed both one and two segments per cell side and 
favored 2 segments per side as opposed to the conventional 1 
segment per side.) The model was run at selected frequen- 
cies from 3 to 450 MHz, in free space. These parameters were 
chosen to span geometry parameters where guidelines for NEC 
wire modeling have been presented. These guidelines were: 


1. Ratio between segment length (D) and segment radius 
(a) greater than or equal to two. 


2. Angle between wires greater than 2m(D/a) “1+ 
3. Spacing between segments greater than or equal to 10a. 
Theoretically, this wire-grid box with the dipole inside 
should radiate. Therefore, average power gain and input 
impedance will be calculated. If acceptable values are 
obtained, a statement concerning the accuracy of near-field 
calculations via NEC can be prepared. In this case, average 
power gain must be close to one (free space). Values of 
acceptable input impedance must be limited to those where 
the ratio of the real to the imaginary part are less than 
104. This limit is because of the accuracy of calculating 
complex numbers in NEC. 
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1 MT. DIPOLE IN 2MT. WIRE GRID BOX 





THETA = 60.00 PHI = 50.00 ETA = 90.00 


Figure 3.1 1 Meter Dipole in 2 Meter Wire- 
Grid Box. 
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C. MONOPOLE ON WIRE GRID BOX 

Bhattacharya, Long and Wilton [Ref. 2] performed an 
experimental investigation to determine the input admittance 
characteristics of a monopole antenna mounted on a con- 
ducting cubical box over a ground plane. In this case, a 6 
cm. monopole was positioned at the center, edge, and corner 
on the top of the box and the admittance was measured (see 
Table 1). This measured data [Ref. 2] was compared with 
calculated data obtained using a five-sided WIRE-GRID BOX 
model of .1 meter per side and cells of .0125 by .0125 


meters. 






Conducting Box 
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Figure 3.2 Geometry of a Monopole Mounted on a Cube Shaped 
Box Over a Ground Plane. [From Ref. 2] 
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TABLE 1 





6 cm. MONOPOLE POSITION 









MEASURED DATA FOR A 6 CM. MONOPOLE ON A CUBE BOX [REF.2] 
- 
GHZ. 
AT EDGE 3.50cm |AT CORNER 5.15cm. 
FROM CENTER FORM CENTER 
|e 


In order to simulate the cube box in Figure 3.2, a model 





AT CENTER 























of a 6 ecm. monopole divided into 5 segments was placed on 
top of the wire-grid box at the center (see Figure 3.3), at 
the edge 3.75 cm. from the center (see Figure 3.4), and at 
the corner 5.3 cm. from the center (see Figure 3.5). The 
feed of the antenna was placed at the base and second 
segment up on the monopole at the center of the wire grid 
only to see if there was any major difference in the 
results. 

The calculated average power gain and input impedance 
are shown in Tables 2, 3, and 4. 
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D. MONOPOLE ON SURFACE PATCH BOX 

In the same fashion as the wire grid box, a five sided 
SURFACE-PATCH BOX of .1 meter per side and cells .02 by .02 
meters was constructed in order to simulate the cubical box 
described in Reference 1. For this model, two different 
tops of the cube were designed to allow placement of the 
model of a 6 cm. monopole divided into five segments at the 
center (see Figure 3.6), at the edge 3.5 cm. from the center 
(see Figure 3.7) and at the corner 4.93 cm. from center (see 
Figure 3.8). In each case, the monopole was fed at the base 
segment. 

The calculated average power gain and input impedance 
are shown in Tables 5, 6, and 7, and description of the 


results are discuss in Chapter V. 
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6 CM. MONOPOLE AT CENTER OF WIRE GRID BOX 





THETA = 60.00 PHI = 60.00 ETA = $0.00 


Figure 3.3 6 cm. Monopole at Center of Wire- 
Grid Box. 
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6 CM. MONOPOLE AT EDGE 3.75 CM. FROM CENTER 





ETA =a58200 


THETA = 60.00 PHI = 60.00 


MOn@pOle wat Hage 6c osseme 


from Center. 
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Figure 3.4 
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6 CM. MONOPOLE AT CORNER 5.3 CM. FROM CENTER 





THETA = 70.00 FHI = 65.00 ETA = 90.00 


Mmucaittens.5 6 Cm. Monopole at Corner 5.3 cm. 
from Center. 
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TABS 22 


CALCULATED DATA FOR A MONOPOLE AT THE CENTER 
OF A WIRE-GRID BOX, BASE SEGMENT FEED 


FREQUENCY |AVG. Power IMPEDANCE / 
Giug. Gain ADMITTANCE 


13.0-j40.7 |15.1-345.4 |14.1-343.0 
7.154+322.2 16.49+519.6 |6.82+320.9 
GHZ. 


16.8-j23.3 |17.9-j25.5 |17.4-j24.4 


20. 3+0/262 2 18.4+}]26.3 19.44927535 


el 4 iGHZeadere 1.93 |21.6-j6.60 |24.1-j7.00 |22.9-j6.80 
42.4+413.0 |38.2+j11.1 |40.3+512.6 


1.15 GHZ.11.82 |1.94 |27.4+479.30 |2922+410.3 |28.3+79me0 
32.7-j11-1 |30.4-5310.7 |31.6—9mloma 
1.2 GHZ.|1.83 |1.96 |34.3+324.1 |36.9+327.4 |35.6+j25.8 
19.5-j13.7 |17.4-j13.0 |18.5-j13.4 
1.3 GHZ./1.84 {1.97 |50.34+350.0 |52.14+356.2 [|51.2+553.1 
10.0-j10.8 |8.88-39.57 |9.44-j10.2 
1.4 GHZ./1.86 |1.98 |66.7+370.6 |66.8+579.0 |66.8+374.8 
7.07-j7.48 |6.24-j7.38 |6.70-j7.43 
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TABLES 3 


CALCULATED DATA FOR A MONOPOLE AT THE EDGE 
OF A WIRE-GRID BOX, BASE SEGMENT FEED 


FREQUENCY |AV. POWER IMPEDANCE/ 
GHZ. GAIN ADMITTANCE 
WO GHZ sie 7 28.9-j45.8 
9.90+515.6 
10S) (GHZ llbee 1.0 31.4-330.8 
eis GHZ 44). 08 34.4-316.0 
Cosa eye 
115 2687 i110 38.1-31.08 
BiG On 4 
i Deere sais 43.0+313.5 
21.2-j6.67 
1.3 GHZ./1.40 55.4+4541.0 
11.7-38.63 
1.4 GHZ.|/1.65 71.1+365.7 
7.60-37.00 
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TABLE 4 


CALCULATED DATA FOR A MONOPOLE AT THE CORNER 
OF A WIRE-GRID BOX, BASE SEGMENT FEED 


FREQUENCY |AVG. POWER | IMPEDANCE/ 
GHZe GAIN ADMITTANCE 
eae GHZ.) 1ess6 48.0- 747.4 
TOM 6) 10.4 

GHZ. 49.9-}34.5 
13.7439 .37 


52.0-j6.63 
16.6+36.61 


1.15. GHZ. laeeey 55.3-36.63 
7c eee 
1.2 GHZ dee 60.1+37.60 
16.4-j2.07 
1.3. GHZ.]1.88 74.3+334.6 
11.1-35.10 
1.4 GHZ.|1.88 92.7+355.0 
8.00-34.73 
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6 CM. MONOPOLE AT CENTER OF SURFACE PATCH BOX 





THETA = 60.00 FHIi= 60.00 EIA = 90.00 


Figure 3.6 6 cm. Monopole at Center of Surface- 
Patch Box. 
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6 CM. MONOPOLE AT EDGE 3.5 CM. FROM CENTER 





THETA = 60.00 FPol= 60 00 aia Seago 


Figure 3.7 6 cm. Monepolesgaemrdgqe s 5ecr- 
from Center. 


oe 


6 CM. MONOPOLE AT CORNER4.93 CM. FROM CENTER 





THETA = 70.00 PHI = 65.CO ETA = $9.00 


Figure 3.8 6 cm. Monopole at Corner 4.93 cm. 
from Center. 
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TABLE 5 


CALCULATED DATA FOR A MONOPOLE AT THE CENTER 
OF A SURFACE-PATCH BOX, BASE SEGMENT FEED 


FREQUENCY |AVG. POWER | IMPEDANCE/ 
GHZ. GAIN ADMITTANCE 
100 “GHZ.4/2..072 1LOs8=s20)..2 

20.6+338.4 


05> “Gaze 14.3-j5.49 
60.6+}323.2 
ia ee aes: 58 
aro = aes: 9 
1.15 GHZai 2c 24.2+321.9 
22.9 4 2.0'.0 
1.175 "GHA 27 2 3 taco 2 
17.1-j18.2 
nee GHZ.|2.01 30.6+9]34.2 
14, 5916.1 
1.225 GHZ eee 34.0+740.0 
1223 =e 
1.309 GHZ 2-08 44.8+)]55.5 
8.84-j310.9 
1.4 GHZ.|2.02 58.4+]71.8 
6.81-]8.37 
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TABLE 6 


CALCULATED DATA FOR A MONOPOLE AT THE EDGE 
OF A SURFACE-PATCH BOX, BASE SEGMENT FEED 


FREQUENCY |AVG. POWER | IMPEDANCE/ 
GHZ. GAIN ADMITTANCE 
GHZ. 22.6+47.40 
39.0-j13.1 


GHZ. 25.3+j18.4 
25.9-j18.8 


ee ILO) (eg, |) IES 2S eo 4 
Wie, Say 1a eS 
17 Oe GH Zemeiee cio 32.9+740.0 
2 a ae 
ile 5 GHZ.) 1). 88 2 re ao © 
Oma — oo 7 
Le GHZ 2) 188 37.0701. > 
O72 5 ale co 
1.225, GHA. | 1.88 40.34357.0 
8277 =) 1 7 
49°.54+972.3 
6.44-]9.41 
Oe ee © 
1.94-}3.46 
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TABEE 7 


CALCULATED DATA FOR A MONOPOLE AT THE CORNER 
OF A SURFACE PATCH BOX, BASE SEGMENT FEED 


FREQUENCY |AV.POWER IMPEDANCE/ 
GHZ. GAIN ADMITTANCE 
1.00 GHZ s> 34.1-]53.2 
8.524713 .3 
1.05 GHZ.|1.86 36.2=940.5 
12.2+j313.7 
1:10, GHZ. |e. s7 38.7-j320.4 
17 Oa) ee eel 
1.15 GHZ |ieses A135 Oege4. 2 
20.9+36.93 


1.175 "GHA. cs £5. Veqieos 
212540. D2 
1.20 GHZ. is? A] -6=ay oo 
2 gO Galea © 
1.225 GHA. eo > 504 aon Ol 
19.5-j2.33 
ies GHZ. |1.83 60.6+524.7 
14 = 75.70 
ee GHZ.|1.82 77.4+345.9 
9.54-j5.66 
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IV. WIRE GRID SOURCE BOX 


The results of the average power gain and input 
impedance of the wire-grid cube-shaped box model are 
displayed in Table 8. For all of the frequencies chosen, 
the average power gain was close to one. (This model was in 
free space where the theoretical average power gain must be 
one). Only input impedance values for frequencies from 42 
to 450 MHz can be considered useable because the ratios of 
the real part to the imaginary part are less than 10%. For 
frequencies of 150 and 300 MHz, where the side dimension of 
the box is one and one-half wavelength, respectively, the 
ratio of the real part to the imaginary part was 109, which 
indicates that the dipole is resonant (A/2 and A long) and 
radiates very well through the widely-spaced wire grids. 
However when the model gets smaller in relation to the 
wavelength, as for 30 MHz and below, the results have ratios 
of 10% and higher between the real and the imaginary parts 
of the input impedance. 

Calculations of the electrical near field for the 
frequencies of 30, 75 and 300 MHz, were performed for a path 
starting at the interior dipole and passing through the wire 
grid to the outside of the box, (see Figures 4.1 and 4.2 and 
Appendix A). The plots show the E field along the Y and Z 


axis. It can be seen in all the plots that the fields 


eo 


inside the box at the wire grid walls are down several 
orders of magnitude from those at the dipole's surface. 
Passing through the grid, the fields drop suddenly by 
another order of magnitude for the two lower frequencies. At 
300 MHz, where the grid spacing is very wide (.4A) the box 


is quite transparent, with little shielding effect. 
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TABLE 8 


CALCULATED DATA FOR THE DIPOLE INSIDE 
THE WIRE-GRID BOX 


FREQUENCY! LAMBDA} AV. POWER IMPEDANCE 
Mhz. Mts. GAIN 
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Figure 4.1 E-Field along Dipole Axis (Z-Axis) for 1M. 
Dipole inside 2M. Wire-Grid Box at 30 MHz. 
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Figure 4.2 E-Field perpendicular to Dipole Axis (Y-Axis) 
for 1M. Dipole inside 2M. Wire Grid Box at 
30 MHZ. 
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V. RESULTS OF THE VARIATION OF THE POSITION 
OF THE MONOPOLE ON THE BOX 

A. WIRE GRID CASE 

To describe the curves which show the measurements and 
the calculation of the input impedance for the variation of 
the position of the 6 cm. monopole, five composite graphs 
were constructed showing the curves of NEC-calculated and 
measured Conductance (G) and Susceptance (B) from Reference 
2, (see Figures 5.1 through 5.5). In the case of the 
monopole at the center of the box, the resonant frequency, 
as defined by the axis crossing of the susceptance curve, 
decreased from 1.138 to 1.125 GHz for Figure 5.1 (base 
segment feed), decreased from 1.138 to 1.122 GHz for Figure 
5.2 (second segment feed), and decreased from 1.138 to 1.123 
GHz for Figure 5.3 (the average of the two different feed 
methods). All these are relative to the measured suscep- 
tance from Reference 2. The variation was small due to the 
feed position. The feed segment was varied only for this 
position of the monopole. In some cases, when the feeding 
segment is tied to many other segments, errors may result. 
When feeding a vertical segment tied to four horizontal 
segments, calculations are difficult because the segments 
which are fed usually like to have a straight segment on 
both sides. Due to the good results obtained with the base 


segment feed, all the other comparisons are with the 


SZ 


monopole base segment feed; the measured data was also 
obtained by feeding that segment. The peak of the real part 
of the input admittance was lower for the calculated value, 
changing from 46.4 to 41 mS. 

In the case of the monopole at the edge of the box, the 
calculated resonant frequency decreased slightly from 1.164 
to 1.155 GHz and the peak conductance varied from 31.7 to 
26.2 mS (Figure 5.4) in comparison to measurements. In the 
case of the monopole at the corner, the NEC resonant 
frequency decreased from 1.192 to 1.157 GHz and the peak 
Gemauctance from 20.9 to 17.8 mS (Figure 5.5). For Sehis 
wire-grid model, the results of the calculated data are very 
close to the measured values from Reference 2. Also, the 
result of the average power gain for this model were close 
to two (correct value of average power gain for perfect 


conducting ground plane). 


B. SURFACE-PATCH CASE 

For this model, the results of the calculated (NEC) 
input admittance were very good for the monopole located at 
the corner. There was a shift in resonant frequency of 
about -5 % for the monopole at center, (see Figure 5.6). For 
the edge-mounted geometry, the correlation between 
measurement and calculation was quite poor (see figure 5.7). 
The calculated values of the average power gain were close 
to two for all positions of the monopole (see Tables 5, 6, 
and 7, Chapter III). For the monopole at the center of the 
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surface-patch box, the resonant frequency decreased from 
1.138 to 1.052 GHz and the peak conductance increased from 
46.4 to 60.6 mS (see Figure 5.6). The resonant frequency 
for corner mounting changed from 1.192 to 1.2 GHz and the 
peak conductance increased from 20.9 to 21.5 mS (see Figure 
5.8). This accuracy might be due to the fact that in this 
case the top of the surface-patch box had been more finely 
divided into smaller patches than for the other two mounting 
cases. This was done in order to position the monopole at a 
specific distance from the center, close to the distance 


used in Reference 2. 
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Figure 5.1 
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Figure 5.2 


36 


(ZHD) AONANDAYS 
Ort Sel O€'l GZ el cit Olt GOt Oot 





Oc- 


Oc- 


OF 





p> 
of 
es 
eae ih = ee aie 
mm 
| a 
oo a 
7 ia 
| (oan) : = 
Seu (Ovals 
GN39D41 . o 
Ba 
© 
$0334 °D45S OGNZ 38 3ASVE JO SJDVHFA 
RS we ad ; 


XO¥ GHD SAYIM 40 YLNOD LY FIOdONOW 


nter of Wire-Grid Box, 
second Segment Feeds. 


6 cm. Monopole at ce 
Average of Base and 


Figure 5.3 


od 


(ZHO) AONSNOSYS 
Ov'l cel oet Zt Oz Gt Olt SOl 


—_—.— — — -——_— 


— (HLND WOH WO G 
{HIN | WOYS WO 


XO GHD SHIM SO ADCS LV AIOdONOW 





(SOHWITHIW) SONVLLINGYV 


6 cm. Monopole at Edge of Wire-Grid 


Box, Base Segment Feed. 


Figure 5.4 
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6 cm. Monopole at Corner of Wire-Grid 


Figure 5.5 
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Figure 5.6 
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VI. CONCLUSION AND RECOMMENDATIONS 


The initial goal of this work was to investigate the 
accuracy of near-field calculations using the Numerical 
Electromagnetic Code (NEC) on ship-like structures. Since no 
validation bench mark results were available for _ ship 
antenna installations, a modeling exercise was undertaken, 
which would yield insight for near-field work in the future. 
Box-like structures were analyzed, one with an interior 
dipole and one with a monopole mounted on the top. For the 
dipole inside the wire grid, there were no measured values 
to compare to, but the calculations of the average power 
gain and input impedance were as expected. The E field 
Variation was also examined as a function of frequency. The 
fields inside the box drop off as the observation point 
moves away from the source. For a non-radiating structure, 
the fields would be similar to the interior of a small 
cavity, with standing wave patterns. Any input resistance 
in the perfect cavity case could only come from dipole 
conductor loss. The leaky nature of the wire box produced 
non-zero input resistance for the dipole, but as_ the 
frequency was lowered to a point where the grid size 
approached »/20, the real part of the input impedance was 
numerically swamped out by the reactance, an expected result 


arising from less radiation leakage. The value of the leaky 
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box model is that it raises concerns about wire grids used 
as models of ships structures for HF RADHAZ predictions. On 
a real-world ship, the antennas/sources will be located on 
the exterior of the ship. The predicted near fields around 
these radiators may be subjected to substantial calculation 
errors if the wire grid surfaces are not "tight" enough. 
Measurements of near-fields for well-controlled generic 
shambened shapes are needed as a validation benchmark for 
additional numerical model developments. 

In the case of the monopole over the cube-shaped box, 
the two modeling methods that NEC has available were tried. 
The wire-grid model worked very well, but the surface patch 
gave good results only for the monopole positioned at the 
corner of the box where the patches at the top are smaller. 
If smaller patches (finer surface current definition) for 
‘the other two positions of the monopole had been used the 
results might have been better. Those data sets are easy to 
produce, but will be difficult to calculate because of 
excessive computer time requirements. It is recommended 
that this be done on a super-computer in the continuation of 
this work. 

This study is an important step in the direction of 
near-field modeling of shipboard shapes because it shows the 
kind of wire-grid definition needed to get accurate input 
impedance for antennas on boxes and gives an idea of how to 


apply surface-patch modeling for accurate input impedance. 
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APPENDIX A 


E-~FIELDS ALONG DIPOLE AXIS (Z-AXIS) AND 


PERPENDICULAR AXIS (Y-AXIS) FOR 1M DIPOLE INSIDE 
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Figure A.1 E-Field Along Dipole Axis (Z-Axis) for 1M. 
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Dipole Inside 2M. Wire Grid Box at 75 MHz. 
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A.2  E-Field Perpendicular to Dipole Axis (Y-Ax1s) 
For 1M. Dipole Inside 2M. Wire Grid Box at 75 MHz. 
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APPENDIX B 
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APPR NDiPx 
DESCRIPTION OF NEC 

The Numerical Electromagnetics Code (NEC) is a user- 
oriented computer code for analysis of the electromagnetic 
response of antennas and other metal structures. It 
built around the numerical solution of integral equations 
for the currents induced on the structure by sources or 
incident fields. This approach avoids many of the simpli- 
fying assumptions required by other solution methods and 
provides a highly accurate and versatile tool for electro- 
magnetic analysis. 

The code combines an integral equation for smooth 
surfaces with one specialized to wires to provide for 
convenient and accurate modeling of a wide range of 
structures. A model may include nonradiating networks and 
transmission lines connecting parts of the structure, 
perfect or imperfect conductors, and lumped element loading. 
A structure may also be modeled over a ground plane that may 
be either a perfect or imperfect conductor. 

The excitation may be either voltage sources on the 
structure or an incident plane wave of linear or elliptic 
polarization. The output may include induced currents and 
charges, near electric or magnetic fields, and radiated 
fields. Hence, the program is suited to either antenna 


analysis or scattering and EMP studies. 
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The integral equation approach is best suited to struc- 
tures with dimensions up to several wavelengths. Although 
there is no theoretical size limit, the numerical solution 
requires a matrix equation of increasing order as_ the 
structure size is increased relative to the wavelength. 
Hence, modeling very large structures may ‘require more 
computer time and file storage than is practical on a 
particular machine. In such cases tanec high-frequency - 
approximations such as geometrical optics, physical optics, 
or geometrical theory of defraction may be more suitable 
than the integral equation approach used in NEC. 

NEC contains the Numerical Green's Function for 
partitioned-matrix solution and a treatment for lossy 
grounds that is accurate for antennas very close to the 
ground surface. It also includes an option to compute 
Maximum coupling between antennas and useful options for 
Emremcteure input.[{Ref.6:p.1] 

1. ZONING CONSIDERATIONS 

NEC is a discrete sampling code where a_ complex 
structure must be discretized into a number of simple 
elements (wires or surface patches) to which the Electrical 
Field Integral Equation (EFIE) or Magnetic Field Integral 
Equation (MFIE) are applied. As with any description of the 
real world, there are approximations, but because of the 


versatility in modeling the geometry of ae structure, 
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its approximations more closely resemble nature. The 
resemblance is strongly influenced by the choice of zoning 
(l.e., dissecting) the structure in the program. The 
smaller the geometric elements, the closer the model comes 
to reality. However, the smaller the elements, the larger 
number of elements, which means the larger the matrix of 
equations and hence, the more costly the solution. There is 
a point beyond which smaller zones do not yield a 
substantially more accurate solution; it may be necessary to 
rerun the problem with increasingly smaller elements to find 
the point of diminishing returns. The choice of proper 
zoning then is gained with experience and becomes as much of 
an art as it is a science. The guidelines for the science 
aspect are as follows. ([Ref.6] 
a. Wires 

Segments should follow the paths of conductors 
using a piece-wise linear fit on curves. Generally, 
segments lengths (f\) should be less than 1X ; shorter 
segments (.05,\) or less may be needed at critical regions 
(junctions or curves). Segments smaller than 1073 ox should 
be avoiced since the similarity of constant and cosine 
components lead to numerical inaccuracy. The radius of the 
wire (a) relative to A depends on the Kernel use in the 
Integral Equation. Two options exist. The thin-wire Kernel 
models a filament current, while the extended Kernel models 


a uniform current distribution around the segment surface. 
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The field of the distributed current is approximated by the 
first two terms in a series expansion of the exact field, in 
powers of a*. The first (a°) term is identical to the thin- 
wire Kernel; the second term extends the accuracy for larger 
values of a. Both Kernels incorporate the thin-wire 
approximations and both require 2ma/,. «1. The thin-wire 
Kernel requires A /a > 8; the extended Kernel relaxes this 
to A/a > 2. These values ensure errors are less than 1%. 
The extended Kernel is used at free wire ends and between 
parallel, connected segments. The thin-wire Kernel is 
always used at bends. 
is Surfaces 

A conducting surface is modeled by small flat 
surface patches which conform as closely as possible to 
curved surfaces. The parameter defining a patch is a normal 
unit vector, originating from the center of the patch, 
defined in Cartesian coordinates. Each patch must be 
connected by a wire at the patch center for the program to 
integrate the surface current. The code divides each patch 
into four equal patches about the wire end, along the unit 
vector lines describing the surface of the of the patch. An 
interpolation function is applied to the four patches to 
represent placement of the current onto the patches, and the 
function is numerically integrated. Patches with wires 
connected to them (active patches) should be chosen to be 


approximately square with sides parallel to the unit vectors 
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defining the surface. Only one wire may connect to a patch, 
that wire may not be connected to another patch. A minimum 
of about 25 patches should be used per wavelength of surface 
area; the maximum size for an individual patch is about .04 
square wavelength. The number of patches used increases, 
and the size of each patch decreases, as the radius of 
curvature decreases. Smaller patches should be used at 
edges since the current amplitude may vary rapidly in this 
region. Long narrow patches should be avoided. Patches are 
restricted to modeling voluminous bodies with closed 
surfaces; parallel surfaces on opposite sides cannot be too 
close together. 
ec: Ground Plane 

For a perfectly conducting ground, the code gen- 
erates a reflected image. Structures may be close to, or 
contact, the ground; however, for a horizontal 
wire: /h? +a? >107© where a = wire radius, h = height of wire 
aX1s above the ground plane, and h/a 2 3. 

A finitely conducting ground may be modeled by an 
image modified by the Fresnel Plane-wave Reflection 
coefficients. This method is fast, but of limited accuracy 
and should not be used for structures close to the ground, 
or having a large horizontal extent over the ground. The 
Sommerfeld/Norton model uses the exact solution and is 
accurate close to the ground; the horizontal restriction is 


the same as for a perfect ground. This model is only used 
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for wire to wire interactions, for surfaces the code reverts 


to Fresnel Reflection coefficients. [Ref. 6 ie 


Do 


APPENDIX D 
RADIATION HAZARD (RADHAZ) ON BOARD SHIPS 
The near fields of Navy antennas are primarily 
important as radiation hazards (RADHAZ). There are essen- 
tially three types of radiation hazards: 
RADHAZ to Personnel (HERP) 
RADHAZ to Ordnance (HERO) 
RADHAZ to Fuel (HERF) 
These hazards and their general near field criteria are 
discussed in the following sections. 


a. Hazards of Electromagnetic Radiation to 
Personnel (HERP) 


Most Navy platforms possess systems and equipment 
which radiate energy in some portion of the electromagnetic 
spectrum. The radiation can constitute a hazard for 
personnel. The degree to which personnel are exposed to the 
radiation depends on the type of platform, the systems 
installed, the location of the systems (particularly their 
antennas), the possibility of directive antennas illumin- 
ating areas occupied by personnel, etc. 

The effects of exposure to this radiation can be 
classified into two categories, direct and indirect. Direct 
radiation is defined as the energy impinging directly on the 
body. Indirect radiation is defined as the capture of EM 
energy by a metallic object. Resulting voltage and currents 
cause burns when the metallic objects are contacted by 
personnel. 
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With certain reservations, physiologists generally 
agree that damage to personnel from direct radiation is 
primarily a heating effect. 

In recent years, however, many distinctly 
nonthermal effects have been documented, some of which have 
been show to be dependent on peak powers whose average value 
is not great enough to produce heating. Frequency 
dependence, with no heating, has also pnareeterized many of 
the observed effects. While the full significance of these 
effects as human hazards has not been established, the fact 
that they occur at average power levels considered to be 
negligible suggest that, at the least, an awareness of their 
existence should be assumed. Some recorded nonthermal 
effects listed in MIL-HDBK-238(NAVY) 1973 are: 

(1) Minor changes in human blood properties upon exposure 
to EM energy of proper frequency and intensity. 

(2) Auditory response. Certain people hear a buzz when 
exposed to microwave radiation. The sensation of 
sound is’ probably not the microwave frequency but 
response to the pulse repetition frequency. 

(3) Abnormalities of the chromosome structure occurring 
upon exposure. 

(4) Movement, orientation, and polarization of protein 


molecules in pulse RF fields. 
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(5) Unexplained response of man to radar. Epigastric 
distress, emotional upsets, and nausea may 
occasionally occur at as low as 5-10 mW/cm? and are 
most commonly associated with the frequency range 
from 

8 X 10° to 12 X 103 MHz 

Changes in the transport rate of materials across the 
blood-brain barrier in humans and animals have also been 
observed. The significance of these changes is not yet 
understood. 

The Navy RADHAZ standard follows the current ANSI 
standard for 3 to 30 MHz with constant values below 3 MHz. 
Personnel shall not be exposed to a power density which, 
when averaged over any 0.1 hour period, exceeds 
900/f? (MHz) mW/cm’? in the frequency domain of 3 to 30 MHz. 
Neither the root mean squared electric field strength (E) 
nor the root mean squared magnetic field strength (H) may 


exceed the following values when averaged over any 0.1 hour 


period: 
1897 
E> -tmnz 7% 
4.74 
H=—- waz 4/™ 


(These are the electric and magnetic field strengths roughly 
corresponding to electromagnetic waves in free space to 
which the value of power density show above may be as- 


Signed.) For a condition where exposure is not regular in 
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time or continuous in level over the 0.1 hour period, the 
equivalent energy fluence level of one-tenth of the power 
density number in units of mW-h/cm? may be use as the limit 
of exposure for any 0.1 hour period. In situations where 
measurements of two or more quantities are available, the 
most restrictive shall be used as the limiting factor. 

RF burns constitute an indirect radiation effect. 
Ship communication transmitters induce voltages on various 
metallic structures such as underway replenishment gear, 
booms, and loading hooks. The voltages are capable of 
Causing painful burns to personnel who come in contact with 
these structures. 

It has been established that 140 volts or greater, 
measured from the object to ground, can cause burns to 
persons touching the object. 


b. Hazards of Electromagnetic Radiation to Ordnance 
(HERO) 


Hazards of electromagnetic radiation to ordnance 
(HERO) stem from the use of sensitive electroexplosive 
devices (EEDs) in ordnance systems. The principal emphasis 
is on protecting electroexplosive devices such as squibs, 
detonators, explosive switches, and ejection cartridges. 
Premature actuation of some of these devices can, of course, 
have dire consequences, such as rocket ignition or warhead 
detonation on a flight deck. Premature actuation by RF can 
also impact the reliability of a system, as for example when 
the power supply for a missile guidance system is expanded 
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inadvertently by the RF initiation of a squib switch. In 
addition, EED firing characteristics can he altered by 
electromagnetic induced heating. 

For the above reasons NAVMAT Instruction 5101.1 
requires that weapons systems and devices containing EEDs be 
reviewed and tested (if deemed necessary) and positive 
certification obtained that they can be handled with 
impunity in the aera predicted electromagnetic environ- 
ment before introduction for service use. Three classifica- 
tions pertinent to HERO for ordnance items have been 
established. The classifications are based on weapons 
susceptibility. The degree of susceptibility is dependent 
upon the electromagnetic environment, the potential for the 
induction of electromagnetic energy into the ordnance 
system, and the characteristic of the EED. Items that are 
negligibly susceptible and require no field intensity 
restrictions beyond general requirements during all phases 
of normal employment are classified as HERO SAFE ORDNANCE. 
Items that are moderately susceptible and require field 
intensity restrictions for at least some phases of employ- 
ment are classified as HERO SUSCEPTIBLE ONDNANCE SYSTEMS. 
Items that are highly susceptible and require severe field 
intensity restrictions for some or all phases of employment 
are classified as HERO UNSAFE ORDNANCE. 

An extensive testing program exists to define HERO 


problem areas and ensure the safety and reliability of 
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ordnance items. Tests are conducted in the maximum RF 
environment to which ordnance items will be exposed during 
the stockpile-to-launch sequence. This implies that HERO 
SAFE ORDNANCE items have been tested to be safe and reliable 
in the criteria environments. 


c. Hazards of Electromagnetic Radiation to Fuels 
(HERF) 


The possibility of accidentally igniting gasoline 
at shore facilities and aboard aircraft carriers and other 
ships handling fuel has been considered when rf-indu::ea 
sparks have been observed. Many years ago the Navy 
conducted extensive tests and investigations of ignition of 
fuel vapors by rf-induced sparks and was found that the risk 
of ignition was small because the following conditions must 
occur simultaneously for ignition to take place: 

(1) Fuels must be heated above their flash points. 
Fuels will not ignite Cees oe is provided in 
certain specific and exact proportions. Fuels will 
also not detonate unless certain exact proportions of 
oxygen and fuel are maintained. 

(2) Sufficient energy must be provided to sustain 
ignition. Although radiant energy is the primary 
source of energy associated with HERF phenomena, no 
mechanism exists by which this energy can interact 
directly with the fuel. Radiated energies are capable 
of causing arcing in susceptive structures. It is 
these arcs which produce ignition of the fuels. 
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An arc iS a voltage breakdown between two fixed 
electrodes. It takes a certain amount of voltage to break 
down the dielectric between the electrodes, but it takes 
energy behind this source of voltage to sustain the arc for 
sufficient duration to ignite fuel mixtures. A spark, on 
the other hand, is formed when two metallic objects in 
electrical contact are separated and there is sufficient 
energy between these electrodes to sustain the spark and 
thus cause the ignition. This is similar to the situation 
in which a fueling nozzle is withdrawn from the tank opening 
when an aircraft is being fueled. The radiation field 
intensity in the proximity of a radiating antenna could be 
sufficient to induce the required energy. 

From actual measurements, it has been determined 
that a spark of energy of 50 V-A is required to ignite 
gasoline in an explosive vapor test device. Recently, some 
attempt has been made to relate the fuel hazard to electric 
field intensity, the primary result has been to show that 
the igniting electric field intensity is a function of 
frequency and iS minimal in the upper hf band. 

In recent years, the probability of occurrence of 
this particular hazard has been further decreased with the 
following advances of technology: 

(1) Ship Design. In the design of ship topside 
arrangements, the location of fuel handling stations, 


fuel vents, etc, 1s considered relative to the 
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(2) 


(3) 


placement of RF transmitting antennas. Every effort 
is made to locate these antennas away from the fuel 
stations and vents. 

Fuel Nozzle Modifications. Heat shrinkable tubing is 
applied to nozzles used in conventional gasoline 
fueling. This prevents metal-to-metal contact between 
the nozzle and the filler tube, thus eliminating the 
possibility that a spark will occur. Another 
hazard-reducing development in the aircraft refueling 
area has been the introduction of the pressurized 
fueling approach. In this procedure, the fuel nozzle 
contacts the metallic filler tube of the aircraft 
prior to the opening of the valve to the tank. Thus, 
if a spark should occur, there would be no flammable 
fuel-air mixture to ignite. 

Fuel. The type of fuel now being most extensively 
used is JP=-5, which has low volatility compared to 
gasoline. Because JP-5 does not have the volatility 
of gasoline, it is not likely to produce flammable 
fuel-air mixtures under accidental circumstances, this 
has done much to reduce the fuel hazard aboard _ ships. 
The fuel oil which powers the ships is still less 


flammable than JP-5.[Ref. 4:pp. 4-9] 
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